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Abstract 

We develop robust pricing and hedging of a weighted variance swap when market prices for a finite 
number of co-maturing put options are given. We assume the given prices do not admit arbitrage and 
deduce no-arbitrage bounds on the weighted variance swap along with super- and sub- replicating 
strategies which enforce them. We find that market quotes for variance swaps are surprisingly close 
to the model-free lower bounds we determine. We solve the problem by transforming it into an 
analogous question for a European option with a convex payoff. The lower bound becomes a problem 
in semi- infinite linear programming which we solve in detail. The upper bound is explicit. 

We work in a model-independent and probability-free setup. In particular we use and extend 
FoUmer's pathwise stochastic calculus. Appropriate notions of arbitrage and admissibility are intro- 
duced. This allows us to establish the usual hedging relation between the variance swap and the 'log 
contract' and similar connections for weighted variance swaps. Our results take form of a FTAP: we 
show that the absence of (weak) arbitrage is equivalent to the existence of a classical model which 
reproduces the observed prices via risk-neutral expectations of discounted payoffs. 

Key Words: Weighted variance swap, weak arbitrage, arbitrage conditions, model-independent 
bounds, pathwise Ito calculus, semi-infinite linear programming, fundamental theorem of asset pric- 
ing, model error. 
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1 Introduction 



In the practice of quantitative finance, the risks of 'model error' are by now universally appreciated. 
Different models, each perfectly calibrated to market prices of a set of liquidly traded instruments, may 
give widely different prices for contracts outside the calibration set. The implied hedging strategies, 
even for contracts within the calibration set, can vary from accurate to useless, depending on how well 
the model captures the sample path behaviour of the hedging instruments. This fundamental prob- 
lem has motivated a large body of research ranging f rom a s set pr i cing through portfoli o optimisation, 
risk m a nagement to macroeconom ics (see for example [Con] (|2006f ) . iFollmer et al.l (|2009t) . lAcciaio et al.l 
(2011), Hansen and SargentI ( 20101 ) and the references therein). Another stream of lite rature, t o whic h 
this paper is a contribut i on, develops a robust a pproach to mathematical finance, see iHobsonI (|l998t ) , 
Davis and HobsonI (120071 ) . ICox and Obloil (|2011al) . In contrast to the classical approach, no probabilistic 
setup is assumed. Instead we suppose we are given current market quotes for the underlying assets and 
some liquidly traded options. We are interested in no-arbitrage bounds on a price of an option implied 
by this market information. Further, we want to understand if, and how, such bounds may be enforced, 
in a model-independent way, through hedging. 

In the present paper we consider robust pricing and hedging of a weighted variance swap when prices 
of a finite number of co-maturing put options are given. If (St, t g [0, T]) denotes the price of a financial 
asset, the weighted realized variance is defined as 



(1.1) 



RVT = Y.h{Su) (log 
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where i,; is a pre-specified sequence of times = to < ti < ■ ■ ■ < tn ~ T , in practice often daily sampling, 
and h is a given weight function. A weighted variance swap is a forward contract in which cash amounts 
equal to A X RVt and A x P^^ are exchanged at time T, where A is the dollar value of one variance 
point and is the variance swap 'price', agreed at time 0. Three representative cases considered in 
this paper are the plain vanilla variance swap h{s) = 1, the corridor variance swap h{s) ~ where / 

is a possi b ly sem i-infinite interval in M+, and the gamma swap in which h{s) = s. The reader can consult 
Gatherall tOOd) for information about variance swaps and more extended discussion of the basic facts 



presented below. In this section we restrict the discussion to the vanilla variance swap; we return to the 
other contracts in Section 

In the classical approach, if we model 5*4 under a risk-neutral measure Q af0 St = Ft exp(X( — ^ {X)t), 
where Pi is the forward price (assumed to be continuous and of bounded variation) and Xt is a continuous 
martingale with quadratic variation process {X)t, then St is a continuous semimartingale and 

log ^ = {g{t^) - + {Xt, - Xt,_J, 

with g{t) = log(Pt. ) — ^(X)t.. For m = 1, 2 ... let {s™, i = 0, . . . , fc,„} be the ordered set of stopping times 
in [0,r] containing the times ti together with times r™ 0,t™ = inf{i > rjp_-^^ : \Xt - Xrjji^ \ > 2"™} 
wherever these are smaller than T. If RV™ denotes the realized variance computed as i n (11.11) but using 



the times s™, then RVp {X)t = {^ogS)T almost surely, see iRogers and Williamsl (|2000f ). Theorem 
IV. 30.1. For this reason, most of the pricing literature on realized variance studies the continuous-time 
limit {X)t rather than the finite sum (jl.l|) and we continue the tradition in this paper. 

The key insig ht into the analysis of variance derivatives in the continuous limit was provided by 
iNeubergeil (1994). We outline the arguments here, see Section 2] for all the details. With St as above and 



/ a function, it follows from the general Ito formula that Yt = .f{St) is a continuous semimartingale 
with d{Y)t = (/'(5i))^d(5)f . In particular, with the above notation 



1 .„ 1 .,„, 1 .„ 1 



d(log50 = -dSt~:^d{S)t^-dSt--d{logS)t, 

so that 



St 252 ' St 2 



(1.2) (log S}t^2 [ ^dSt - 2 log(5T/^o) 

Jo "^t 

This shows that the realized variation is replicated by a portfolio consisting of a self-financing trading 
strategy that invests a constant $2Dt in the underlying asset S together with a European option whose 
exercise value at time T is X{St) = —2log{ST/ Ft). Here Dt is the time-T discount factor. Assuming 
that the stochastic integral is a martingale, (jl.2[) shows that the risk-neutral value of the variance swap 
rate P^^ is 

(1.3) P^" = E[(log5)Tl = -2IE[log(5T/PT)] = 2-^Piog, 

that is, the variance swap rate is equal to the forward value Piog/P'r of two 'log contracts' — European 
options with exercise value — log(S'T/PT), a convex function of St- (|1.3p gives us a way to evaluate P^^ 
in any given model. The next step is to rephrase Piog in terms of call and put prices only. 

Recall that for a convex function / : R+ M. the recipe f"{a,b] = f^{b) — f^{a) defines a positive 
measure f"{dx) on £?(M+), equal to f"{x)dx if / is . We then have the Taylor formula 

I^Xf, poo 

(1.4) fix) = fixo) + r+{xo)ix~xo)+ {y-x)+v{dy)+ (x - y)+ f" {dy) . 

Jo JxQ 

Applying this formula with x ^ St,xq = Ft and /(s) = log(s/5o), and combining with (jl.2p gives 

(1.5) l{\ogS)T= r ^dSt-log{FT/So) + l-ST/FT+ ^^^—^dK+ H i^l^-pldK. 



2' ^ Jo St ' ..V . ' • Jo Jpr 

Assuming that puts and calls are available for all strikes K S R"'" at traded prices Pk, Ck, (|1.5p provides 
a perfect hedge for the realized variance in terms of self-financing dynamic trading in the underlying (the 



-"^For example, St = {e^^ ''''5o)(e'^^* 2" *) in the Black-Scholcs model, with the conventional notation. 
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first four terms) and a static portfolio of puts and calls and hence uniquely specifies the variance swap 
rate P^^ as 

2 f°° 1 

(1-6) ^D^l ]^(^^^(^^^-) + CKliK>Fr))^K. 

Our objective in this paper is to investigate the situation where, as in reality, puts and calls are 
available at only a finite number of strikes. In this case we cannot expect to get a unique value for P^^ 
as in (|1.6p . However, from (|1.2[) we expect to obtain arbitrage bounds on from bounds on the value 
of the log contract — a European option whose exercise value is the convex function —log{ST/FT)- It 
will turn out that t he weighted varianc e swap s we consider are associated in a similar way with other 



convex functions. In iDavis and HobsonI (j2007h conditions were stated under which a given set of prices 



(Pi,...,P„) for put options with strikes {Ki, . . . , Kn) all maturing at the same time T is consistent 
with absence of arbitrage. Thus the first essential question we have to answer is: given a set of put 
prices consistent with absence of arbitrage and a convex function At , what is the range of prices Pa for 
a European option with exercise value At(5't) such that this option does not introduce arbitrage when 
traded in addition to the existing puts? We ans wer this question in Sect ion [3] following a statement of 
the standing assumptions and a summary of the IPavis and HobsonI (j2007t) results in Section [2j 

The basic technique in Section [3] is to apply the duality theory of semi-infinite linear programming 
to obtain the values of the most expensive sub-rcplicating portfolio and the cheapest super-replicating 
portfolio, where the portfolios in question contain static positions in the put options and dynamic trading 
in the underlying asset. The results for the lower and upper bounds are given in Propositions 13.11 and 
13.51 respectively. (Often, in particular in the case of plain vanilla variance swaps, the upper bound 
is infinite.) For completeness, we state and prove the fundamental Karlin-Isii duality theorem in an 
appendix. Appendix [Xj With these results in hand, the arbitrage conditions are stated and proved in 
Theorem 13.61 

In this first part of the paper — Sections [5] and [5] — we are concerned exclusively with European options 
whose exercise value depends only on the asset price St at the common exercise time T. In this case 
the sample path properties of {St,t £ [0,T]} play no role and the only relevant feature of a 'model' 
is the marginal distribution of St- For this reason we make no assumptions about the sample paths. 
When we come to the second part of the paper. Sections U and [SJ analysing weighted variance swaps, 
then of course the sample path properties are of fundamental importance. Our minimal assumption is 
that {St,t G [0,T]} is continuous in t. However, this is not enough to make sense out of the problem. 
The connection between the variance swap and the 'log contract', given at (jl.2p . is based on stochastic 
calculus, assuming that St is a continuous semimartingale. In our case the starting point is simply a set 
of prices. No model is provided, but we have to define what we mean by the continuous-time limit of the 
realized va riance and by continuous time trading. An appropriate framework is the pathwise stochastic 



calculus of lFollmerl (jl98lh where the quadratic variation and an Ito-like integral are defined for a certain 
subset Q of the space of continuous functions C[0,r]. Then (jl.2p holds for S{-) G Q and we have the 
connection we sought between variance swap and log contract. For weighted variance swaps the same 
connection exists, replacing '— log' by some other convex function A^, as long as the latter is C^. In the 
case of the corridor variance swap, however, AJ^ is discontinuous and we need to restrict ourselves to a 
smaller class of paths C for which the Ito formula is valid for functions in the Sobolev space W2 ■ AH of this 
is consistent with models in which St is a continuous semimartingale, since then ¥{uj : S{-,uj) G C} = 1. 
We discuss these matters in a separate appendix, Appendix [BJ 

With these preliminaries in hand we formulate and prove, in Section 2] the main result of the paper, 
Theorem 14. 3[ giving conditions under which the quoted price of a weighted variance swap is consistent 
with the prices of existing calls and puts in the market, assuming the latter are arbitrage- free. When the 
conditions fail there is a weak arbitrage, and we identify strategies that realize it. 

In mathematical finance it is generally the case that option bounds based on super- and sub-replication 
are of limited value because the gap is too wide. However, here we know that as the number of put options 
increases then in the limit there is only one price, namely, for a vanilla variance swap, the number P^^ 
of (|1.6p , so we may expect that in a market containing a reasonable number of liquidly-traded puts the 
bounds will be tight. In Section [5] we present data from the S&P500 index options market which shows 
that our computed lower bounds arc surprisingly close to the quoted variance swap prices. 

This paper leaves many related questions unanswered. In particular we would like to know what 
happens when there is more than one exercise time T, and how the results are affected if we allow 
jumps in the price path. The latter question is considered in the small time-to-expiry limit in the recent 



paper iKeller-Ressel and Muhle-Karbd ()2011l) . Exact pricing formulas for a wide variety of contracts on 
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discre tely-sampled realized variance, in a general Levy process setting, are provided in Crosby and DavT^ 
(|2nilh . 



2 Problem formulation 

Let {St,t G [0,r]) be the price of a traded financial asset, which is assumed non-negative: St £ M+ = 
[0, oo). In addition to this underlying asset, various derivative securities as detailed below are also traded. 
All of these derivatives are European contracts maturing at the same tim e T. The present time is 0. We 



make the following standing assumptions as in lDavis and HobsonI (120071 ): 



(i) The market is frictionless: assets can be traded in arbitrary amounts, short or long, without transaction 
costs, and the interest rates for borrowing and lending are the same. 

(ii) There is no interest rate volatility. We denote by Dt the market discount factor for time t, i.e. the 
price at time of a zero-coupon bond maturing at t. 

(iii) There is a uniquely-defined forward price Ft for delivery of one unit of the asset at time t. This will 
be the case if the asset pays no dividends or if, for example, it has a deterministic dividend yield. We let 
Tt denote the number of shares that will be owned at time t if dividend income is re-invested in shares; 
then Fq = So and Ft = So/iDtTt) for t > 0. 

We suppose that n put c iptions with strikes < < ... < A'„ and maturity T arc traded, at 
time-0 prices Pi, . . . , P„. In (iDavis and HobsonL[2007l) . the arbitrage relations among these prices were 



investigated. The facts are as follows. A static portfolio X is simply a linear combination of the traded 
assets, with weights tti, . . . , 7r„, 0, ?/) on the options, on the underlying asset and on cash respectively, it 
being assumed that dividend income is re-invested in shares. The value of the portfolio at maturity is 

n 

(2.1) Xt = ^ TT, [K, -St]+ + (J^TtSt + yjD^^ 

i=l 

and the set-up cost at time is 

n 

(2.2) ^0 = ^ ^^P^ + (I^So + V'- 

1=1 

Note that it does not make any sense a priori to speak about Xt - the value of the portfolio at any 
intermediate time - as this is not determined by the market input and we do not assume the options 
are quoted in the market at intermediate dates t e (0,r). A model is a filtered probability space 
(ri,F = {J^t)t£lo,T]:^) together with a positive F-adaptcd process {St)te[o,T] such that 5*0 coincides with 
the given time-0 asset price. Given market prices of options, we say that is a market model for 
these options if Mt = St /Ft is an F-martingale (in particular, St is integrable) and the market prices 
equal to the P-expcctations of discounted payoffs. In particular, is a market model for put options 
if Pi = DT&[Ki — St\^ for i = 1, • • • , n. We simply say that is a market model if the set of market 
options with given prices is clear from the context. It follows that in a market model we have joint 
dynamics on [0,r] of all assets' prices such that the initial prices agree with the market quotes and the 
discounted prices are martingales. B y the (easy part of the) First Fundamental Theorem of Asset Pricing 



(jPelbaen and Schachermaveiv Il994f ) the dynamics do not admit arbitrage. 

Our main interest in the existence of a market model and we want to characterise it in terms of the 
given market quoted prices. This requires notions of arbitrage in absence of a model. We say that there 
is model-independent arbitrage if one can form a portfolio with a negative setup cost and a non-negative 
payoff. There is weak arbitrage if for any model there is a static portfolio (tti, . . . , 7r„, </), "0) such that 
Xq < but ¥{Xt > 0) = 1 and F \Xt > 0] > 0- In par ticular, a model independent arbitrage is a 



special case of weak arbitrage, as in ICox and Obloi (2011b). Throughout the paper, we say that prices 



are consistent with absence of arbitrage when they do not admit a weak arbitrage. 

It is convenient at this point to move to normalized units. If for t = 1, . . . , n we define 

Pi Ki 
(2.3) p^ - 75"^, h - 

then, in a market model, we have Pi ~ DTFT^[Ki/ Ft — Mt]^ and hence 

p^=E[h-MT] + . 
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Figure 1: Normalized put prices {ki,pi) consistent with absence of arbitrage. An additional put with 
price r — and strike k £ [0, a] does not introduce arbitrage. 



Also, we can harmlessly introduce another put option with strike fco = and price po — 0. Let 
n = max{i : pi = 0}, n = inf{^ : p,; = kj — 1 }(= + oo if {••■} = 0) and K = [n, h]. 

The main result of iDavis and HobsonI (|2007l . Thm. 3.1), which we now recall, says that there is no 



weak arbitrage if and only if there exists a market model. The conditions are illustrated in Figure [T] 

Proposition 2.1 Let r : [0, fe„] — ^ be the linear interpolant of the points {ki,pi),i = 0, . . . ,n. Then 
there exists a market model if and only if r[-) is a non-negative, convex, increasing function such that 
r(0) = 0,r(fc) > [fc — 1]+ and r'_(fcfi/\„) < 1, where r'_ denotes the left-hand derivative. If these con- 
ditions hold except that h ~ oo and r'_(fc„) = 1 then there is weak arbitrage. Otherwise there is a 
model-independent arbitrage. Furthermore, if a market model exists then one may choose it so that the 
distribution of St has finite support. 



Remarks, (i) In iDavis and HobsonI (j2007f ). the case of call options was studied. The result stated 
here follows by put-call parity, valid in view of our frictionless markets assumption. The normalised call 
price is Ci — Pi -\- 1 — ki. 

(ii) Of course, no put option would be quoted at zero price, so in applications n ~ always. As will 
be seen below, it is useful for analysis to include the artificial case n>0. 

Throughout the rest of the paper, we assume that the put option prices (Pi, . . . ,Pn) do not admit a 
weak arbitrage and hence there exists a market model consistent with the put prices. Let be a market 
model and let p be the distribution of Mt in this model. Then p satisfies 



(2.4a) 
(2.4b) 
(2.4c) 



/ 1 n{dx) = 1 
x p{dx) = 1 



[h - x]'^p{dx) = Pj, i = 1, 



Conversely, given a probability measure p on which satisfies the above we can construct a market 
model A4 such that p is the distribution of A fy. Fo r exam ple, let (il, F, P, {Wt)t^R+) be the Wiener space. 
By the Skorokhod embedding theorem (cf. (Obloj, 20041 )). there is a stopping time r such that WV ~ fi 



A(t/(T-t)) 



for 



and (WrAt) is a uniformly intcgrable martingale. It follows that we can put Mt = 1 + 
t G [0,T). This argument shows that the search for a market model reduces to a search for a measure p 
satisfying (|2.4p . We will denote by Mp the set of measures fi satisfying the conditions 

Lemma 2.2 For any p e Mp, p{R+\K) = 0. 



Proof. That /i[0, fc„) = when n > follows from (j2.4c[) with i ~ n. When n < n we have Cn = 0, i.e. 
there is a free call option with strike kn and wc conclude that fi{kn, oo) = 0. □ 
The question we wish to address is whether, when prices of additional options are quoted, consistency 
with absence of arbitrage is maintained. As discussed in Section [1] we start by considering the case where 
one extra option is included, a European option maturing at T with convex payoff. 
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3 Hedging convex payoffs 

Suppose that, in addition to the n put options, a European option is offered at price P\ at time 0, with 
exercise value Xt{St) at time T, where At is a convex function. We can obtain lower and upper bounds 
on the price of At by constructing sub-replicating and super-replicating static portfolios in the other 
traded assets. These bounds are given in Sections 13.11 and 13.21 respectively and are combined in Section 
l3.3l to obtain the arbitrage conditions on the price P\- 

We work in normalised units throughout, that is, the static portfolios have time-T values that are 
linear combinations of cash, underlying AIt and option exercise values [ki — Mt]'^ ■ The prices of units 
of these components at time are Dt, Dt and DxPi respectively, where a unit of cash is $1. Indeed, to 
price Mt observe that $1 invested in the underlying at time yields TtSt/Sq = St/FtDt ~ Mt/Dt 
at time T. To achieve a consistent normalization for At we define the convex function A as 

(3.1) \{x) = ^Xt{Ftx). 

In a market model M. we have P\ — Dt^[\t{St)] = DtFt^[\{Mt)], so the normalized price is 

Px = - E[A(A/t)] 

and the cost for delivering a payoff A(A/t) is DtPx- 

3.1 Lower bound 

A sub-replicating portfolio is a static portfolio formed at time such that its value at time T is majorized 
by A(A/t) for all values of Mt- Obviously, a necessary condition for absence of model independent 
arbitrage is that DtPx be not less than the set-up cost of any sub-replicating portfolio. It turns out 
that the options ki with i < n or i > n arc redundant, so the assets in the portfolio arc indexed by 
fc = 1 , . . . , m where 

m = (nH-l)An — n + 1 
and the time-T values of these assets, as functions of x ~ Mt are 

ai(x) = 1 (Cash) 

(3.2) a2{x) = X (Underlying) 

ai+2{x) = [fc„+j - + , 1 = 1,..., TO -2 (Options). 

We let a(x) be the m-vector with components ak{x). Note that a,„(a;) is equal to — x]~^ if n < n 

and to [kn — x]'^ otherwise. The set-up costs for the components in p.2p . as observed above, are Dt, Dt 
and DTPn+i respectively. The corresponding forward prices are as in (|2.4p : 

h = 1 

(3.3) 62 = 1 

bi+2 = Pn+i, i=l,...,m-2. 

We let b denote the m-vector of the forward prices. A static portfolio is defined by a vector y whose kth 
component is the number of units of the kth asset in the portfolio. The forward set-up cost is y'^b and 
the value at T is y^a(AfT)- 

With this notation, the problem of determining the most expensive sub-replicating portfolio is equiv- 
alent to solving the (primal) semi-infinite linear program 

Plb : sup y^b subject to y^a{x) < X{x) \/x G K. 



The constraints are enforced only for x G K. If n > [n < n] we have a free put with strike fc„ [call with 
strike kn] and, since A is convex, we can extend the sub-replicating portfolio to all of K"*" at n o cost. 
The key res ult here is the basic duality theorem of semi-infinite linear programming, due to llsiil 



and Karlin, see Karlin and Studden (Il966l) . This theorem, stated as Theorem lA.il and its proof are given 



in Appendix [Xj The dual program corresponding to Plb is 

I?LB : inf / X(x)ii(dx) subject to / SL(x)fi(dx)—h° 
msmJk Jk 
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where M is the set of Borel measures such that each is integrable. The constraiuts iu Dlb can be 
expressed as /j, satisfying (j2.4p ioi n < i < n. This is simply equivalent to /i G Mp since, as shown in 
Lemma[n21 any G Mp has support in K. Let Vp and be the values of the primal and dual problems 
respectively. It is a general and easily proved fact that Vp < V^. The 'duality gap' is Vj^ — Vp . The 
Karlin-Isii theorem gives conditions under which there is no duality gap and we have existence in Plb- 

Proposition 3.1 We suppose as above that X{x) is a convex function on M^, finite for all a; > 0, and 
that {ki,pi) is a set of normalised put option strike and price pairs which do not admit a weak arbitrage. 
If \{x) is unbounded as x ^ Q and n = then we further assume thatpi/ki < p2/fe- Then V^ = Vp and 
there exists a maximising vector y. The most expensive sub-replicating portfolio of a European option with 
payoff \t {St) at maturity T is the static portfolio as in (j2.1[) with weights i/j^ = FtDtYi, 4>^ = y^/^T, 
TT^+i = y2+i for i = 1, . . . , m — 2 and 7r| = otherwise. For this portfolio, = DtFtV^ . 

If there is existence in the dual problem Dlb then there is an optimal measure /i^ which is a finite 
linear combination of Dirac measures /i^ = X^JLi 'Wj5xj{Ax) such that each interval [kj^ fcj+i) contains at 
most one point Xj . For this measure 

n 

(3.4) /i^({x}) > ^ y tt\[K, - Ftx\+ + (I)^VtFtx + ijj'^ = \t{Ftx). 

St — Ft'J^ — ^ 
i—1 

Proof. The first part of the proposition is an application of Theorem lA.ll The primal problem Plb 
is feasible because any support line corresponds to a portfolio (containing no options). The functions 
ai, . . . ,am are linearly independent. Recall from Proposition 12.11 that if {(fci,pi)} do not admit weak 
arbitrage then there is a measure /i satisfying the conditions (|2.4p and such that is a finite weighted sum 
of Dirac measures. It follows that |A(x)|/z(dx) < oo unless one of the Dirac measures is placed at x = 
and A is unbounded at zero. If n > then there is no mass on the interval [0, n), hence none at zero. When 
n = we always have pi/fci < ^2/^2- If Pi/^i — ^2/^2 then the payoff [fc2 — Mt]'^ — ^2(^1 — Mt\^ /ki has 
nuU cost and is strictly positive on (0, ^2). Since pi > there must be some mass to the left of fci, and this 
mass must be placed at 0, else there is an arbitrage opportunity. But then J = +cxd and V^ = +cxd. 
The condition in the proposition excludes this case. In every other case there is a realizing measure /i such 
that ^({0}) = 0. Indeed, if pi/fci < ^2/^2 then the extended set of put prices {(fc, 0), (fci,pi), . . . , (fc„,p„)} 
is consistent with absence of arbitrage if fc £ [0, a], where a = (fcip2 — k2Pi)/iP2 —pi) (see Figure [1}. Any 
model realizing these prices puts weight on the interval [0, k). Thus V^ is finite under the conditions we 
have stated. It remains to verify that the vector b belongs to the interior of the moment cone M,„ defined 
at (|A.ip . For this, it suffices to note that for all i such that ki £ (fc„, kn) it holds that [ki — 1]+ < pi < ki, 
and so the condition is satisfied. We now conclude from Theorem lA.ll that Vp = V^ and that we have 
existence in the primal problem. The expressions for ^'^ etc. follow from the relationships (j2.3l) between 
normalized and un-normalizcd prices. 

Assume now that the dual problem has a solution. Any optimal measure £ Mp satisfies 

X(x)^^(dx) = inf < / \(x)ij,(dx) 

Recall K = [fc„, kn] and partition K into intervals . . . , InAn+i defined by 

h = [h-i,ki) for i = n + 1, . . . 7l A n and InAn+i = [knAn, kn), 

so IfiAn+i = if ?i < Lemma 13.21 below asserts that we may take atomic with at most one atom in 
each of the intervals Ii. By definition the optimal subhcdging portfolio X^ satisfies X^, < Xt{St) while 
our duality result shows that the expectations of these random variables coincide. Hence X^ = Xt{St) 
a.s. for Mt ~ St /Ft distributed according to ij,^ , and p.ip follows. □ 

Lemma 3.2 Let fi G Mp and suppose J-^ \X{x)\fj,{dx) < 00. Define 

(3.5) If, ^ {i <n An + l\jj.{h) > 0}. 
Now let ji' be the measure 



in which Ox denotes the Dirac measure at x, and for an index i Q I, Xi — — '^^^^ — . Then fi G 



and 



(3.6) / X{x)fi'{dx) < / X{x)n{dx). 
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Proof of lemma: The inequality p.6p follows from the conditional Jensen inequality. A direct com- 
putation shows that n' satisfies (|2.4p . so that ^' G Mp. □ 

It remains now to understand when there is existence in the dual problem. We exclude the case when 
At is affine on some [z^oo) which is tedious. We characterise the existence of a dual minimiser in terms 
of properties of the solution to the primal problem and also present a set of sufficient conditions. Of the 
conditions given, (i) would never be encountered in practice (it implies the existence of free call options) 
and (iii),(iv) depend only on the function Xt and not on the put prices Pi. The examples presented in 
Section |3. 1 . II show that if these conditions fail there may still be existence, but this will now depend on 
the Pi. Condition (ii) is closer to being necessary and sufficient, but is not stated in terms of the basic 
data of the problem. 

Proposition 3.3 Assume Xt is not affine on some half-line [z,oo). Then, in the setup of Proposition 
the existence of a minimiser in the dual problem Dlb fails if and only if 

(3.7) n = oo and = (jiTTS + ip"^ Dj^^ < Xt{s), for all s > Kn- 

St — s 

In particular, each of the following is a sufficient condition for existence of a minimiser in -Dlb ■ 
(i) n < oo; 
(ii) we have 

(3.8) X| = (P^TtK,, + i^^D^r,^ < XriKn) and lim Xt{s) - fli^Trs = oo; 

ST=Kn S-S-OO 

(Hi) for any y < XxiKn) there is some x > Kn such that the point {Kn,y) lies on a support line to Xt 
at x; 

(iv) Xt satisfies 

/•oo 

(3.9) / xX'Tidx) = +00 . 

Jo 

Proof: We consider two cases. 

Case 1: fi < n, i.e. condition (i) holds. In this case the support of any measure fi G Mp is contained 
in the finite union In+i U ■ ■ ■ U of bounded intervals, and X]r=ri+i l^i-^i) = 1- Further, Pn-i > kn — 1 
and (|2.4b[) together imply that /i(/fi) > 0. Let j = 1, 2, ... be a sequence of measures such that 



X{x)fj,j{dx) ^ inf < / X{x)^{dx)> as j — > oo. 

By Lemma |3.2[ we may and do assume that each fXj is atomic with at most one atom per interval. We 
denote = A*j(-^i) and let Xj denote the location of the atom in Ii. For definiteness, let Xj = A, where 
A is some isolated point, if fij has no atom in Ii, i.e., if it;* ~ 0. Let A be the set of indices i such that 
converges to A as j — ?> oo, i.e. ^ A for only finitely many j, and let B be the complementary set of 
indices. Then there exists j* such that J2i£B ~ ^ ^'^^ ^ > J* ■ Because the and Xj are contained in 
compact intervals there exists a subsequence jk, k = 1,2... and points w\, x\ such that w*^ — >■ w\ and 
x\ xl, as k ^ oo. It is clear that y^, .-nwl = 1 and that 

Vt^\\ui f X{x)fi,,{dx) = f X{x)fiHdx) 

where fj.^{dx) ~ J2ieB ^l^xii^^)- Similarly, the integrals of l,x and [ki — x]+ converge, so fi^ G Mp. 
Finally, since the intervals Ii are open on the right, it is possible that xl G li+i. If also x^"*"^ G li+i we 
can invoke Lemma 13.21 to conclude that this 2-point distribution in li^i can be replaced by a 1-point 
distribution without increasing the integral. Thus retains the property of being an atomic measure 
with at most one atom per interval. We have existence in the dual problem and, by the arguments above, 
xl = Xt{St) dn^-a.s. 
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Case 2: n = oo. Let X^{s) denote the payoff of the portfoho X'^ at time T when 5*^ = s. Suppose that 
a minimiser exists in Dlb, which we may take as atomic. Since n = oo there exists x € (fcn, cxd) such 
that fJ,^{x}) > and hence, by ((5^ . xI,{Ftx) = X{Ftx). This shows that ([XT)) fails since Ftx > K^. 

It remains to show the converse: that if n = oo but (j3.7p fails then a minimiser fi^ exists. Given 
our assumption on At, if p.7p fails then :— sup{s > : x\,{s) — At(s)} € [Kn, oo). We continue 
with analysis similar to Case 1 above. Here we have intervals In, ■ ■ ■ , In covering K„ = [kn, kn), but also 
a further interval In+i = [kn, oo) which does not necessarily have zero mass. Since is unbounded, 
the argument above needs some modification. We take a minimizing sequence of point mass measures 
Hj as in Case 1. If x^^^ converge (on a subsequence) to a finite x""*"^ then we can restrict our attention 
to a compact [0,x"'^^ + 1] and everything works as in Case 1 above. Suppose to the contrary that 
liminf x"^^ = oo, as j — > cxo. We first apply the arguments in Case 1 to the sequence fij of restrictions 
of 11 j to K„. Everything is the same as above except that now Wj < 1. A subsequence converges to a 
sub-probability measure /i^ on K„, equal to a weighted sum of Dirac measures as above. Define fi^ = /i^ 
if Li = /it(K„) = 1 and otherwise /i^ = /i^ + (1 — ii)Sx, where x S In+i is to be determined. Whatever 
the value of x, satisfies conditions (|2.4ap and (|2.4cp (the put values depend only on /i^'). 

Since fij e Mp, J2i=i ^]^) = ^ particular 

n 

l-Y.w]x]>nf;^'kn. 

i=l 

Taking the limit along the subsequence we conclude that 

(3.10) 1 -J^wlxl >knil-Y,<] = ^"(1 - 'l)- 

i=l \ 1=1 / 

By the convergence argument of Case 1, i2 = /o^" 3^/i^(da::) < 1. If ii < 1 we have only to choose 
x = {1 — i2)/(l — ti) to ensure that the 'forward' condition (|2.4b[) is also satisfied. The inequality p.lOp 
guarantees that x > kn- Thus fi^ G Mp and, since n — oo, x > kn and hence > Kn- 

We now show that the complementary case ii = 1 contradicts K'' 6 [i^„,oo). Indeed, if ti = 1 we 
have, since n = oo, 

kn - I <Pn = / [kn - x\^ p) {dx) = / (fc„ - x)p\Ax) = " V] W^X^ . 

Jo Jo 

Observe that then 

n n 

(3.11) = 1 - E '^^^ - — ' 1 - E « = p„ - fc„ + 1 > 0. 



1— J-oo 

i=l i=l 



Since G [/v„, oo), taking 7 A^(A"t + 1) - x|'(A'^ + 1) = A^(X^ + 1) - c/j^Tt > 0, we have 

(3.12) Xt{s)-XI{s)>j{s-K^ + 1), Vs>/v^ + 1. 

Define now a new function At by 

At(s) = AT(s)l,<Kt + (0trTS + 4'^D^^)l,yK^ 

so that we have X'' < Xt{St) < Xt{St) and, by definition, for St > A'^ we have X'' — Xt{St) < Xt{St)- 
It follows that X'' is also the most expensive subreplicating portfolio for At(s) and hence the primal and 
dual problems for A and for X{x) = -^XriFTx) have all the same value. Writing this explicitly and using 
(|3T2)) and (|3ll|) gives: 

0= lim / (A(x)-A(a;))/^,(d.T)= lim(A(x"+i)-A(x"+i)X+i 

> lim 7(.t"+i - - l)w"+i = 7(p„ - fc„ + 1) > 0, 

j — yoo 



a contradiction. 
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We turn to showing that each of (i)-(ro) is sufficient for existence of a dual minimiser /if. Obviously 
n < oo is sufficient as observed above. We now show that either of (ii) or (iii) implies that ift g [Kn, oo) 
and that p.l2p holds. In the light of the arguments above this will be sufficient. X^{s) is linear on 
[Kn,oo) and At is convex hence the difference of the two either converges to a constant or diverges 
to infinity. In the latter case the difference grows quicker than a linear function, more precisely p.l2p 
holds. The former case is explicitly excluded in (ii) and is contradictory with (iii) as the line £{s) := 
(p^Txs + tp'^D^^ + limu^oo {^t{u) — X^{u)) is asymptotically tangential to Xt{s) as s — ?> oo and hence 
the condition in (iii) is violated for any y < i{Kn)- In particular, it follows that < Kn then we can 
add to a positive number of call options with strike Kn and obtain a subreplicating portfolio with 
an initial cost strictly larger (recall that n = oo) than X'^ which contradicts the optimality of X'' . This 
completes the proof that either (ii) or (iii) is sufficient for existence. 

Finally, we argue that (iii) and (iv) arc equivalent. In (iii) the point x satisfies 

Xt{x) + ^{Kn — x) = y, for some ^ G [\'rp{x—),\'r]^{x+)\. 

This equation has a solution x for all y < XxiKn) if and only if lim2;_j.oo Xt{x) — xX'rp{x) = — oo. The 
equivalence with (iv) follows integrating by parts 

/ xX'^idx) = (xX'Tix) ~ Xt{x)) °° . 

□ 

3.1.1 Examples v^rith one put option 

We consider now examples with just one put option price is specified. We illustrate different cases 
when existence in the dual problem holds or fails. For simplicity assume all prices arc normalised, i.e. 
Dt = Ft = 1, and we are given only a single put option with strike k = 1.2. The convex function takes 
the form A(a;) = l/x + ax^ . Computation of the most expensive sub- hedging portfolio can be done by a 
simple search procedure. 

Consider first the case a = 0, so that A(a;) = 1/x. The results are shown in Figure [2] with data shown 
in Tabled) p is the put price, xq,xi the points of tangency, wq^wi the implied probability weights on 
Xq, xi and ^, </> the units of, respectively, cash and forward in the portfolio. 

p = 0.4 is a 'regular' case: we have tangent lines at xq^xi and the solution to the dual problem puts 
weights 8/9, 1/9 respectively on these points. As p increases it is advantageous to include more puts in 
the portfolio, so xi increases. At p = 0.6 we reach a boundary case where "0 = </> = 0, and the put is 
correctly priced by the Dirac measure with weight 1 at = 0.6 = fc/2. Obviously, this measure does 
not correctly price the forward, but it docs correctly price the portfolio since (j) ~ Q. When p > 0.6, the 
only way to increase the put component further is to take ip < (and then clearly the optimal value of 
(j) is 0.) When p = 0.7 the optimal value is tp = —0.8 and we find that in this and every other such case 
the implied weight is wq ^ 1, as the general theory predicts. 



p 


Xo 


Xi 


value 







Wo 


Wl 


0.4 


0.75 


3 


1.2222 


0.6667 


-0.1111 


0.8889 


0.1111 


0.6 


0.6 




1.6667 








1.00 




0.7 


0.5 




2.00 


-0.8 





1.00 





Table 1: Data for Figure [D 



Next, take a = 0.25, 6 = 1, so that A(x) = l/a;-|-0.25x, has its minimum aX x = 2 and is asymptotically 
linear. We take p = 0.7. Taking the expectation with respect to the limiting measure gives the value 
of the lower sub-hedge in Figure [3J but this is not optimal: we can add a maximum number 0.25 of call 
options, which have positive value, giving the upper sub- hedge in Figure [3) This is optimal, but does not 
correspond to any dual measure. 

Finally, let a = 0.0625, 6 = 2, giving A(a;) = l/x + 0.0625x^. The minimum is still at x = 2 but A 
is asymptotically quadratic. This function satisfies condition (iii) of Proposition 13.31 and there is dual 
existence for every arbitrage- free value of p. The optimal portfolio for p = 0.7 is shown in Figure 21 



^In fact for this part of the Proposition we do not need to impose any additional conditions on Ay apart from convexity. 
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Figure 2: Most expensive subhedging portfolios for X{x) = 1/x given a single put option with strike 1.2 
and prices p = 0.4, 0.6, 0.7 




Figure 3: The most expensive subhedging portfolio (upper dashed line) for X{x) ~ 1/x + 0.25.T given 
a single put option with strike 1.2 and price p ~ 0.7. The lower dashed line is the portfolio priced at 
/ X{x)p''{dx) which is suboptinial. 



3.2 Upper bound 

To compute the cheapest super-replicating portfolio we have to solve the linear program 

Pub : inf y^b subject to y^a(a;) > X{x) \/x G K. 

yGR™ 

The corresponding dual program is 

Dub: sup / X{x)iJ,{dx) subject to / a{x)ij{dx) — h, 
tj.eM Jk Jk. 

where, by Lemma we may replace K by M+. By (|1.4p we have for x G M"*" 

A(a;) = A(l) + A'(l+)(a;-l)+ / [k ^ x] + X" (dk) + [ [x ~ k] + X" (dk) . 

•J(0,1] "'(l,oo) 

Consider p e Mp and let p^{k) ~ j[k — a:]^/i(dx), c^(fc) = j[x — k]^ p{dk) be the (normalised) prices of 
puts and calls. Integrating the above against gives 



(3.13) 



X{x)p{dx) = X{l)+ p^ik)X"{dk)+ c^{k)X"idk). 

J(0,1] "'(1,00) 
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Recall that c^(A;) = p^(fc) + 1 — fc and hence maximising in c^(A:) or in Pfj,{k) is the same. p^,{k) is a convex 
function dominated by the linear interpolation of points {ki,pi), i = 0, 1, . . . ,n, extended to the right of 
(kn,Pn) with slope 1. Sin ce we assume the given put prices do not admit weak arbitrage, it follows from 
( Davis and Hobsoill2007t) (see also Proposition l2 . ll above) that this upper bound is attained either exactly 
or in the limit. More precisely, if n = oo one can take pz G Mp supported on {kn, fc„, z}, for z large 
enough, which attain the upper bound on [0,fc„] and asymptotically induce the upper bound on (fc„,oo) 
as z — !• oo. It follows from p.l3p that the value of the dual problem is Vj^ ~ liniz-i-oo A(a:;)/X2(dx). 
If n < n one can take z ~ kn and fik~ attains the upper bound. It follows from (|3.13|) that then 
VE ^ J^XixhkAdx).^ 

From this observations, one expects that y - the solution to the primal problem - will correspond to 
a (normalised) portfolio y''a(a;) which linearly interpolates (fc,;, A(fc,;)), i = fc„, . . . ,n An and (if n = oo) 
extends linearly to the right as to dominate A(fc). The function x t-^ y'^a.{x) is piecewise linear with a 
finite number of pieces and no such function can majorize the convex function A over R"'" unless A(0) < oo 
and A'(oo) = 7 < oo. In general we impose: 

(a) n > or [n — and A(0) < oo) and 
(3.14) ) ' 

(b) n < n or I n = oo and A'(oo) < oo 1 . 

We have the following result. 

Proposition 3.4 // condition (|3.14p holds then there exists a solution y to the linear program P\jb- The 
function y'^a(a;) is the linear interpolation of the points {kn, X{kn)), , . . . , {knAn, ^iknAn)) together with, 
if n = oo, the line l{x) = y^a(fc„) + (x — /c„)A'(oo) for x > kn- Primal and dual problem have the same 
value Vjj ~ Vjj' and the existence of a maximiser in the dual problem fails if and only if n — oo and A 
is not affine on [kn, oo). 

// the condition (j3.14p is not satisfied, there is no feasible solution and = oo. 

Proof. As argued above, (|3.14p is a necessary condition for existence of a feasible solution. Suppose, 
for example, that n = oo = A'(oo). Then (z — kn)fJ.z{{z}) = c„ and hence zfiz{{z}) — c„ as z — > oo. 
Together with A'(oo) ~ oo this implies that Vjj' — oo. Other cases are similar. 

Suppose (|3.14p holds and first consider the case when n = oo. A is bounded on [kn, kn] and the linear 
interpolation is well defined as is the extension beyond Further, there exists some constant 6 such 
that y^a.{x) — X{x) < 6 for all x G R+. The weight y-i+j on the jth put option is the change in slope 
at kn+j, the 'underlying weight y2 is equal to ^, and at x = kn we have A^ 

= yi + y2kn, so the 'cash' 

weight is yi = A„ — fc„7. The value of the objective function is, by definition, 

y^b = / y'^a.{x)p{dx). 
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for any fi E Mp. In particular, since y''a(fci) = X{ki), taking z large enough, we have 
= / y-a(x)Ai.(dx) = 5] A(fe,)/i.({fc.}) + y"a(zK({z}) 

X{x)fi,{dx) + (y"a(z) - A(z))m.({4) < / H^)^^z{dx) + 6pi{{z}). 



(3.15) 



Recall that J x^z{dx) ~ 1 and in particular /i({z}) — > as z — > cxd. Taking the limit as z — >■ oo in the 
above, we conclude that y^b" < . The basic inequality < Vp between the primal and dual values 
then implies that y is optimal for Pub and = Vp . The existence of the solution to -Dub fails unless 
there exists z > kn with y'^a(z) = A(z), which happens if and only if A(z) is affine on [fc„, oo). 

When n <n the arguments are analogous, except that now we need to ensure y'^a(a;) > A(x) only for 
a; G K = [fc„, kn]. The dual problem has a maximiser as observed in the remarks above the Proposition. 
The primal problem also has a solution y but it is not unique. Indeed, let y° = {—km 1, 0, . . . , 0, 1) and 
observe that y^'^b = and y''^a(x-) = for a; G K. In consequence, we can add to y multiples of y° 
without affecting its performance for Pub- D 

We can now summarize the results for the cheapest super-replicating portfolio as in (|2.1[) . (|2.2p . The 
difference with the above is that we need to ensure super-replication for all possible values of St and not 
only for St € [Km Kfi]. Since the payoff Xt, as a function of St, is piecewise linear with a finite number 
of pieces it is necessary that At satisfies 

(3.16) At(0) < oo and \'t{<xi) = 7 < oo. 

Under this condition, the above results show that the cheapest super-replicating portfolio has a payoff 
which linearly interpolates {Ki, XT{Ki)), i = 0, . . . ,n and extends to the right of Kn with slope A^(oo). 

Proposition 3.5 // p.l6p holds then there is a cheapest super-replicating portfolio (i/) *,</)*, tt*) of the 
European option with payoff Xt {St) whose initial price is 



X^ = supDt XT{FTx)fi^{dx). 

2 JR+ 



The underlying component is (jf — ^Tt, the cash component is ip* = DT{XT{Kn) — 7A"„) and the option 
components are 

^ Xt{K,+,) - Xt{K,) _ Xt{K,) - Xt{K,-i) 
If condition (|3.16p is not satisfied, there is no super-replicating portfolio. 



3.3 Arbitrage conditions 

With the above results in hand we can state the arbitrage relationships when a European option whose 
exercise value at T is a convex function Xt{St) can be traded at time at price P\ in a market where 
there already exist traded put options, whose prices Pi are in themselves consistent with absence of 
arbitrage. Recalling the notation of Propositions 13 . 1 1 and 13.51 Xq and Xq are respectively the setup costs 
of the most expensive sub-replicating and cheapest super-replicating portfolios, with Xq = -\-oo when no 
super-rephcating portfolio exists. 

Theorem 3.6 Assume the put prices do not admit a weak arbitrage. Consider a convex function Xt and 
suppose that if Xt is affine on some half-line [z, 00) then it is strictly convex on [0, z). The following are 
equivalent: 

1. The prices Pa, Pi, . . . , P„ do not admit a weak arbitrage. 

2. There exists a market model for put options in which P\ = Dt'^[Xt{St)\, 

3. The following condition (|3.17p holds and either P\ G {Xq,Xq), or P\ = Xq and existence holds in 
DhB, or P\ = Xq < 00 and existence holds in I?uB- 
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K2 

(3.17) P2 > -f-^Pi if Ik — ^ o,i^d X is unbounded at the origin. 

If ((5Tfl) holds and Px i [Xq , Xjf ] then there is a model-independent arbitrage. If (|3.17p holds and 
Px = or Xq and existence fails in -Dlb or in Dtjb respectively, or if p.l7p fails, then there is a weak 
arbitrage. 

Remark 3.7 We note that the robust pricing and hedging problem solved above is essentially invariant 
if At is modified by an affine factor. More precisely, if we consider a European option with payoff 
X^{St) = Xt{St) + <PSt + ^ and let P^i denote its price then the prices Pi, . . . ,Pn,Px o-fc consistent 
with absence of arbitrage if and only if Pi, ... , P„, Px^ = P\ + + ^^Dt are. 

Proof. Suppose first that condition p.l7[) holds. We saw in the proof of Proposition l3.1l tliat tliis condition 
(under its equivalent form p2 > {k2/ki)pi) guarantees the existence of a sub- replicating portfolio with 
value X^. If Pa £ (-'^oi-'^o) then there exists e > such that Px G {X^ + e,XQ - e) and, since 
there is no duality gap, there are measures ^11,^^2 G Mp such that Dt^^-^^[\t{St)] < X^ + e and 
DtE^s [^t{St)\ > — e. a convex combination fi of /ii and ^2 then satisfies Dt^^[Xt{St)] = P\ and 
one constructs a market model, for example by using Skorokhod embedding as explained in Section [5] 
above. If existence holds in Dlb then it was shown in the proof of Proposition 13.11 that the minimizing 
measure fi^ satisfies 

xI = Dt I XTiFTx)ii\dx), 

so that if Px = Xq then jj,^ is a martingale measure that consistently prices the convex payoff At and the 
given set of put options. The same argument applies on the upper bound side. 

Next, suppose that condition p.l7p holds and Px = but no minimizing measure /i exists in I?lb- 
Let A4. he a model and fi the distribution of St under A4. We can, at zero initial cost, buy At(5't) and 
sell the portfolio X^. and this strategy realizes an arbitrage under M if fl{{XT{ST) > X^,}) > 0. Suppose 
now that //({At('S't) > X^}) = and consider two cases. First, if At is not affine on some [z,oo) then, 
by Proposition 13. 3[ fi — 00 and Xt{St) > for St > Kn so that in particular fi[[Kn,oo)) = 0. A 
strategy of going short a call option with strike A"„ (which has strictly positive price since n = cxd) gives 
an arbitrage since St < Kn a-s. in Ai. Second, suppose that X'^{x) = for x > 2: but X'^{x) > for x < z. 
If fi{[Kn, 00)) = then we construct the arbitrage as previously so suppose this does not hold. Recall the 
atomic measure /i^ defined in Case 2 in the proof of Proposition 13.31 and that ti = 1 as wc do not have 
existence of a minimiser for the dual problem. X{x) in p.ip is strictly convex on [0, z) and linear on [z, 00) 
with z = z/ Ft. It is not hard to see that /i^ has to have an atom in some x* £ (fc„_i, fc„) and that either 
z = x* or else z > k„ . Otherwise we could modify fi^ to obtain a minimiser for the dual problem. Strict 
convexity of At on [0, z) implies that there exist at most n points si, . . . , s,„ such that Si £ {Ki-i, Ki) 
and Xt{St) strictly dominates X^ for other values of St < z and hence /i([0,2:)) = /t({si,...,s„}). It 
follows that support of /i^ is a subset of {si/Ft, ■ . . , s„/Ft = .t* }. Let ip > and consider a portfolio Y 
with 

n 

Yt = Yt{St) = ^^'^^'^^ ^ ^t)"^ + ^^'^^ Yt{sj) =0, i = l,...n. 

4=1 

This uniquely specifics 7^ G M. The payoff of Y is simply a zigzag line with kinks in Ki, zero in each Si 
and equal to if) for St > Kn- It follows that (i{Yt > 0} = 1 and /^{It > 0} > as /i([A'„,oo)) > 0. 
However /i^ prices all put options correctly so that the initial price of Y is 



>o = ^ 7»^» + Dt^^ = ^ 7* / {K^ - FTm)^i^ (dm) + Dt^ 

=Dt '^^^^'^ ~ Prmf + j Ai^(dm) ^ YT{FTm)ii\dm) = 0, 



by construction of Y since jj) {{si/ Ft, . ■ . , Sn/FT}) = 1 as remarked above, and where we used ti = 1. 
It follows that Y is an arbitrage strategy in A4. 

Now suppose that condition p.l7p holds and Px = Xq < 00 and there is no maximising measure in the 
dual problem I?ub- Then, by Proposition 13.41 fi = 00 and Xt{St) < X^ for St > Kn. Straightforward 
arguments as in the first case above show that there is a weak arbitrage. 
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Finally, if Pa < ^ model independent arbitrage is given by buying the European option with payoff 
Xt{St) and selling the subheding portfolio. This initial cost is negative while the payoff, since 
subhedges Xt{St); is non- negative. If P\ > Xq we go short in the option and long in the supcrhcdgc. 

Now suppose (|3.17p does not hold, so that P2/K2 ~ Pi/Ki (this is the only case other than p.l7|) 
consistent with absence of arbitrage among the put options). Consider portfolios with exercise values 

H^{St) = [K2 - St]+ - -^[Ki - St]+ 

H2{St) = Xt{St) - \t{K2) - \o{K2){St -K2)-^ (Pa - A(X2) - KiFr - K2)) [K^ - 5t] + , 

where Aq denotes the left derivative. The setup cost for each of these is zero, and Hi{s) > for s G (0, K2) 
while -^2(5) ^ 00 as s 0. There is a number 9 >0 such that H{s) = 0Hi{s)+H2{s) > for s G (0, A'2). 
Weak arbitrage is realized by a portfolio whose exercise value depends on a given model M and is specified 
via 

X (s ) = / "[^'2 - St]+ if P[5t G [0, K2)] = 
'^^ ^' \h{St) if ^[St e [0,K2)] > 0. 

This completes the proof. □ 



4 Weighted variance swaps 



We come now to the second part of the paper where we consider weighted variance swaps. The main 
idea, as indicated in the Introduction, is to show that a weighted variance swap contract is equivalent 
to a European option with a convex payoff and hence their prices have to be equal. The equivalence 
here means that the difference of the two derivatives may be replicated through trading in a model- 
independent way. In order to formalise this we need to define (continuous) trading in absence of a model, 
i.e. in absence of a fixed probability space. Naturally, we face technical difficulties as we need to define 
pathwise stochastic integrals. 

One possibility is to define stochastic integrals as limits of discrete sums. The resul ting object may de- 
pend o n the sequence of partitions used to define the limit. This approach was used in lBick and Willinger 
who interpreted the difference resulting from different sequences of partitions as broker's method 
of implementing continuous time trading order. They were then interested in what happens if they ap- 
ply the pricing-through-replication arguments on the set of paths with a fixed (cr^) realised quadratic 
va riation and wanted t o reco ver Black-Scholes pricing and hedging. However for our purposes the ideas 
of Bick and Willinger (Il994h are not suitable. We are interested in a much wider set of paths and then 
the replication of a weighted variance swap using tradi ng and a Euro pean option would depend on the 



'broker' (i.e. seq uence of partiti ons used). Instead, as in iLvonj (|l995h . we propose an approach inspired 
by the work of iFollmeii (|l98lf ) . We restrict the attention to paths which admit quadratic variation 
or pathwise local time. For such paths we can develop pathwise stochastic calculus including Ito and 
Tanaka formulae. As this subject is self-contained and of independent interest we isolate it in Appendix 
IbI Insightful discussions of this topic are found in iBick and Willinger (Il994h and lLvonj (|l995h . 



To the standing assumptions (i)-(iii) of Section [2] we add another one: 

(iv) {St '■ t < T) E - the set of strictly positive, continuous functions on [0,r] which admit a finite, 
non-zero, quadratic variation and a pathwise local time, as formally defined in Definitions IB.llB.3l and 
Proposition IB .41 of Appendix iBl 

Thus, our idea for the framework, as opposed to fixing a specific model A^, is to assume we arc given 
a set of possible paths for the price process: {St : t < T) E V. This could be, for example, the space 
of continuous non-negative functions, the space of functions with finite non-zero quadratic variation or 
the space of continuous functions with a constant fixed realised volatility. The choice of V is supposed 
to reflect our beliefs about characteristics of price dynamics as well as modelling assumptions we are 
willing to take. Our choice above, P = £+, is primarily dictated by the necessity to develop a pathwise 
stochastic calculus. It w ould b e inter esting to understand if an appropriate notion of no-arbitrage implies 
(iv). A recent paper of I Vovkl (|201lh . based on a game-theoretic approach to probability, suggests one 
may exclude paths with infinite quadratic variation through a no-arbitrage-like restriction, an interesting 
avenue for further investigation. 
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We introduce now a continuous time analogue of the weighted reahzed variance (jl.ip . Namely, we 
consider a market in which, in addition to finite family of put options as above, a w-weighted variance 
swap is traded. It is specified by its payoff at maturity T: 

(4.1) RV^ - P^"(") r w{St/Ft)d{\og S)t - 

Jo 

where Pp^*-^-* is the swap rate, and has null entry cost at time 0. Our assumption (iv) and Proposition lB.6l 
imply that {log St, t <T) € C. Theorem IB . 51 implies that (|4.ip is well defined as long as w € L^^^, we can 
integrate with respect to St or Mt and obtain an Ito formula. This leads to the following representation. 

Lemma 4.1 Let w : — !■ [0, oo) be a locally square integrable function and consider a convex 
function A„, with A", (a) = The extended ltd formula (|B.1[) then holds and reads 

(4.2) A„(Mt) = A^(l)+ / X'^{Mu)dMu + ]- I w{Mu)d{\nM)^. 

Jo ^ J{O.T\ 

The function A^, is specified up to an addition of an affine component which does not affect pricing 
or hedging problems for a European option with payoff \w , see Remark 13.71 above. In what follows we 
assume that w and A„, arc fixed. Three motivating choices of as discussed in the Introduction, and 
the corresponding functions A^,, are: 

1. Realised variance swap: w = \ and Atu(.T) = — ln(x). 

2. Corridor variance swap: w{x) = l(o.a)(2;) or w{x) = l(a.oo)(3;): where < a < co and 

\w{x) = {- In (^) + ^ - l) w(x). 

3. Gamma swap: w[x) = x and A^(a;) ~ x\vi(x) — x. 

Clearly, (|4.2p suggests that we should consider portfolios which trade dynamically and this will allow 
us to link w- weighted realised variance RVj^ with a European option with a convex payoff Xw Note 
however that it is sufficient to allow only for relatively simple dynamic trading where the holdings in the 
asset only depend on asset's current price. More precisely, we extend the definition of portfolio X from 
static portfolios as in (|2.ip - (|2.2p to a class of dynamic portfolios. We still have a static position in traded 
options. These are options with given market prices at time zero and include n put options but could 
also include another European option, a weighted variance swap or other options. At time t we also hold 
Tt4>{Mt) assets St and ipt/Dt in cash. The portfolio is self-financing on (0,r] so that 

(4.3) := <f>iMo)So + ^(0, So) + f 0(Af„)dM„ - Tt<f>iMt)StDt, t e (0, T], 

JO 

and where (j) is implicitly assumed continuous and with a locally square integrable weak derivative so 
that the integral above is well defined, cf. Theorem IB. 51 We further assume that there exist: a linear 
combination of options traded at time zero with total payoff Z = Z{St '■t<T),& convex function G and 
constants (f>, ^ such that 

(4.4) Tt^{Mt)St + iJtlDt>Z-G{Mt)/Dt + ^TtSt+^/Dt, yt<T. 

Such a portfolio X is called admissible. Observe that, in absence of a model, the usual integrability of 
Z is replaced by having finite price at time zero. In the classical setting, the admissibility of a trading 
strategy may depend on the model. Here admissibility of a strategy X may depend on which options are 
assumed to trade in the market. The presence of the term G{Mt) on the RHS will become clear from the 
proof of Theorem 14.31 below. It allows us to enlarge the space of admissible portfolios for which Lemma 
Sil] below holds. 

The two notions of arbitrage introduced in Section [2] are consequently extended by allowing not only 
static portfolios but possibly dynamic admissible portfolios as above. All the previous results remain 
valid with the extended notions of arbitrage. Indeed, if given prices admit no dynamic weak arbitrage 
then in particular they admit no static weak arbitrage. And for the reverse, we have the following general 
result. 



Arbitrage bounds for prices of options on realised variance 



17 



Lemma 4.2 Suppose that we are given prices for a finite family of co-maturing option^. If a market 
model M. exists for these options then any admissible strategy X satisfies ¥\DtXt\ < ^o- In particular, 
the prices do not admit a weak arbitrage. 

Proof. Let 7W be a market model and X be an admissible strategy. We have Xt = +Yt, where 
is a linear combination of payoffs of traded options and Yt = Tt4'{Mt)St + ipt/Dt satisfies (|4.4p . Using 
it follows that 

DtYt = (/)(l)5'o + ^'(0, So) + f 0(M„)dM„ > DtZ - G{Mt) + 0+5'oMt + ^. 

Jo 

We may assume that G > 0, it suffices to replace G by G+. is a market model and in particular Z 
is an intcgrablc random variable. Since the traditional stochastic integral and our pathwisc stochastic 
integral coincide a.s. in Ai, see Theorem IB. 71 we conclude that DtYt is a local martingale and so is 
Nt := DtYt — (j)'^SoMt. We will argue that this implies ENt < Nq. Let p„ be the localising sequence 
for N so that ENt/-,p^ = A^o- In what follows all the limits are taken as n ^ oo. Fatou's lemma shows 
that EN^ < lim inf EA^j"J^p^ and KN^ < lim inf EiVj^p^ . By Jensen's inequality the process G{Mt) is a 
submartingale, in particular the expectation is increasing and EG(M() < 'E,G{Mt) which is finite since 
is a market model. Using the Fatou lemma we have 

limEG(A/fAp„) < EG(Aft) Eliminf G(MtAp„) < liminf EG(AftAp„), 

showing that limEG(MtAp^) = EG(7\/t). Observe that Nt is bounded below by Z — G{Mt), where Z is 
an integrable random variable and G is convex. Using Fatou lemma again we can write 

EG(M()+ + ^Z- - liminf EA^,:;,^^^ = liminf E [G(MtAp„)+ + Z- - iVjAp J > E [G(A/t)+ + - Nt" 

which combined with the above gives EiVj~ = lim inf EA^^^^^ and in consequence EA^t < lim inf EA^tAp„ = 
A^o, as required. This shows that E[DtYt] < Yq. Since in a market model expectations of discounted 
payoffs of the traded options coincide with their initial prices it follows that E[DtXt] < Xq < 0. In 
particular if Xt > and Xq < then Xt = a.s. and the prices do not admit a weak arbitrage. 

□ 

Having extended the notions of (admissible) trading strategy and arbitrage, we can now state the 
main theorem concerning robust pricing of weighted variance swaps. It is essentially a consequence of 
the hedging relation (j4.2p and the results of Section [3] 

Theorem 4.3 Suppose in the market which satisfies assumptions (i)-(iv) the following are traded at time 
zero: n put options with prices Pi, a w -weighted variance swap with payoff (j4.1[) and a European option 
with payoff Xt (St) — FtXw{Mt) and price P\^ . Assuming the put prices do not admit a weak arbitrage, 
the following are equivalent 

1. The 'option' prices (European options and weighted variance swap) do not admit a weak arbitrage. 

2. Pi, . . . , P„, P\,„ do not admit a weak arbitrage and 



(4.5) p-W = ^„2A„,(l). 



3. A market model for all n + 2 options exists. 

Remark 4.4 It is true that under (|4.5p a market model for Pi, . . . ,P„,Py'^'™'' exists if and only if a 
market model for Pi, . . . , P„, Px^ exists. By Theorem \3.6l this is yet equivalent to Pi, . . . , P„, P\^ being 
consistent with absence of arbitrage. However it is not clear if this is equivalent to Pi, . . . , Pn, Pt''^^ 
being consistent with absence of arbitrage. This is because the the portfolio of the variance swap and 
dynamic trading necessary to synthesise —\t{St) payoff may not be admissible when Xt{St) is not a 
traded option. 



These could include European as well as exotic options. 
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Proof. We first show that 2 =^ 3. Suppose that Pi,... ,P„,Pa„ do not admit a weak arbitrage and 
let be a market model which prices correctly the n puts and the additional European option with 
payoff Xt{St). Note that, from the proof of Theorem l3.6[ M exists and may be taken to satisfy (i)-(iv). 
Consider the Ito formula (g^l) evaluated at r„ A T = inf{i > : Mt ^ {l/n,n)} A T instead of T. The 
continuous function is bounded on (l/ri, n), the stochastic integral is a true martingale and taking 
expectations we obtain 



IE[A,„(Af,„AT)] = A^(l) + iE 



w{Mu)d{M)u 



Subject to adding an afSne function to X^, we may assume that > 0. Jensen's inequality shows 
that Xw{Mr„/\T), n > 2, is a submartingale. Together with the Fatou lemma this shows that the LHS 
converges to K[Xyj(MT)] aa n —¥ oo. Applying the monotone convergence theorem to the RHS we obtain 

E[A^(A/t)] = A„,(l) + [RV^] , 
where cither both quantities are finite or infinite. Since is a market model for puts and Xt{St), we 



= 



have E[A^(A/t)] = P\^/{DtFt). Combining (|43)) with the above it follows that E 
and hence is a market model for puts, Xt{St) and the w-wighted variance swap. 

Lemma (j4.2[) implies that 3 1. We note also that, if we have a market model J\A for all the puts, 
Xw{St) option and the w-wighted variance swap then by the above (|4.5p holds. Then Lemma (|4.2p also 
implies that 3 2. 

It remains to argue that 1 2 i.e. that if Pi, . . . , P„, Px^ are consistent with absence of arbitrage 
but (j4.5p fails then there is a weak arbitrage. Consider a portfolio X with no put options, -0(0,50) — 
2Dt{X^{1) - 5oA^(l)) + DtPt"'-'"^ and 0(to) = 2DTX'Jm). The setup cost of X is Xq = Dt{Pt^'''"'^ + 
2A„,(1)). By assumption, A^, is with A^(.t) = w{x)/x'^ 6 L^^^. so that can apply Theorem IB. 51 Then, 
using (|4.3p and (|4.2p we obtain 



Xt - Tt<j)iMt)St + ^t/Dt - 0(5o)S'o/A + 0(0, So)/Dt + ^ [ 0(M„)dM„ 

M Jo 



_ Dt 

~ 15 

= ^A.(M,) - ^ / u,(5„)d(lnM)„ + ^P^^^-^ . 

Observing that 1 > Dt/Di > Dt and /jp w{Su)d{\n M)u is increasing in < it follows that both X and 
—X are admissible. Suppose first that 

(4.6) P-W^ 2^_2A.(1). 

Consider the following portfolio Y: short 2 /Ft options with payoff At(S't); long portfolio X and long a 
w-weighted variance swap. Y is admissible, the initial cost is 

Yo = -2Pa„/Pt + ^0 = -2Px^/Ft + Dt{Pt"'^'"'^ + 2A^(1)) < 0, 

while It = and hence we have a model independent arbitrage. If a reverse inequality holds in (j4.6p 
then the arbitrage is attained by —Y. 

□ 



5 Computation and comparison with market data 
5.1 Solving the lower bound dual problem 



When one of the sufficient conditions given in Proposition l3.3l is satisfied, for existence in the lower bound 
dual problem, then this problem can be solved by a dynamic programming algorithm. We briefiy outline 
this here, referring the reader to Raval ( 201Cl( ) for complete details. For simplicity, we restrict attention 
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to the practically relevant case n = 0, fi = oo. Of course, once the dual problem is solved, the maximal 
sub-hedging portfolio is immediately determined. 
The measure to be determined satisfies 



\(x)Li^(dx) — inf 



and we recall from Lemma 13.21 that we can restrict our search to measures of the form fi(dx) = 

w'i<^xi(d2;) where Xi e [fc^-i, h), Wi > 0, J2i "^i = 1- We denote Co = and, for i > 1, = Si 
the cumulative weight on the interval [0, fc^). For consistency of the put prices ri, . . . , r„ with absence of 
arbitrage. Proposition 12. II dictates that 



and 



for 1 < z < 



,1 



Given C = (Ci, . . . , Cn) (the final weight is of course Wn+i ~ I ^ Cn), the positions Xi are determined by 
pricing the put options. We find that when Ci-i < Ci 

Xi = Xi(Cj-i>Ci) = fcj H ^ ^ lor i^l ,n, 



Xn+l = Xn+liCn) = K + 



Q - Q-1 

1 + r,i fc„ 

1-Cn 



The measure corresponding to policy C is thus 

n 

(5-1) E(C» - C^-l)Sx.iC,-uC,) + (1 - Cn)5x.+l(C„)- 

i=l 

It follows that the minimisation problem inf^gMp Jj^ X{x)fi{dx) has the same value as 
(5.2) 



inf 



C^ti I " C^-l)HX^{C^-lX^)) + (1 - Cn) A(x„+1 (C„)) 



We can solve this by backwards recursion as follows. Define 

K(Cn) = (l-Cn)Xn+l(Cn) 



(5.3) 



VjiCj) = info+iGA,+i,o+i>o{(0+i - Ci)A(x»(Ci, Ci+i)) + Vj+iiQ+i)}, j = 7i - 1, . . . , 0. 



Then Vb(0) ~ v^. For a practical implementation one has only to discretize the sets Aj, and then (|5.3p 
reduces to a discrete-time, discrete-state dynamic program in which the minimization at each step is just 
a search over a finite number of points. 



5.2 Market data 

The vanilla variance swap is actively traded in the over-the-counter (OTC) markets. Wc have collected 
variance swap and European option data on the S&P 500 index from the recent past. Using put option 
prices, the lower arbitrage-bound for the variance swap rate is computed in each case, and summarised 
in Table [5] One sees that the traded price of the variance swap frequently lies very close to the lower 
bound. Nonetheless, under our standing assumptions of frictionless markets, all but one of the prices 
were consistent with absence of arbitrage. We mention that the crash in October 2008 of the S&P 500, 
and indeed the financial markets in general, gave rise to significant increase in expected variance, which 
can be seen in the cross-section of data studied. One data point, the 3-month contract on 20/12/2008, 
lies below our lower bound and appears to represent an arbitrage opportunity. Practitioners have told us 
that this was a day of extreme disruption in the marke10. It is a positive point of our method that we 
are able to pick up such periods of market dislocation, just from the raw price data. 

^Indeed, the final column of Table [2] showing the number of traded put prices used in the calculations, shows that 
December 19 and 20 were far from typical days. 
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Table 2: Historical variance swap (VS) quotes for the S&P 500 index and the lower bound (LB) for 
it, implied by the bid prices of liquid European put options with the same maturity. The units for 
the variance price and LB are volatility percentage points, 100 x yJP^ , and M stands for months. The 
European option price data is courtesy of UBS Investment Bank, and the variance-swap data was provided 
by Peter Carr. 



Term 


Quote date 


VS quote 


LB 


# puts 


2M 


21/04/2008 


21.24 


20.10 


50 


2M 


21/07/2008 


22.98 


22.51 


50 


2M 


20/10/2008 


48.78 


46.58 


93 


2M 


20/01/2009 


52.88 


47.68 


82 


3M 


31/03/2008 


25.87 


23.59 


42 


3M 


20/06/2008 


22.99 


21.21 


46 


3M 


19/09/2008 


26.78 


25.68 


67 


3M 


19/12/2008 


45.93 


45.38 


112 


3M 


20/12/2008 


45.93 


65.81 


137 


6M 


24/03/2008 


25.81 


25.34 


33 


6M 


20/06/2008 


23.38 


23.20 


38 



A The Karlin & Isii Theorem 

Let ai, . . . , arm f be real- valued, continuous functions on K C and let M denote the collection of all 
finite Borel measures on K fulfilling the integrability conditions Jj^ |ai(x)|/^(dx) < oo for 1 = 1,...,™. 
For a fixed vector b = . . . , and letting a(x) = (ai(x), . . . , am(x))'^ for x G K. consider the 

optimisation problem 

(P) : sup y^b s.t. y^a(x) < /(x) Vx e K. 

yeR™ 



Now, define 



(D) : inf / /(x)M(dx) s.t. / a(x)A*(dx) = b. 



where the constraint should be interpreted as /j^ ai(x)/x(dx) = bi, i ~ l,...,m. The values of the 
problems (P) and {D) will respectively be denoted by V{P) and V{D). Finally, Mm C M™ will denote 
the moment cone defined by 

(A.l) l^h= (bi,...,bmV\k^ J a,(x)^(dx), i = l,...,m, A* e m| . 

Theorem A.l Suppose 

1. ai, . . . , a„i are linearly independent over K, 

2. b is an interior point of Mm and 

3. V{D) is finite. 

Then V{P) — V{D) and (P) has a solution. 

Remarks (i) The beauty in the Karlin & Isii theorem is that the pro of draws upon n o more than a 



HEMARKS [1) ine beauty m tnc J\arnn &; isu tncorcm is tnat tne pro oi araws upon n o more tnan a 
finite dir aensional separating hypcrp lane theorem. Proofs can be found in iGlashofj |l979t) and Karlin & 
Studden Karlin and Studde 3 (Il966i . Chapter XII, Section 2). We follow the latter below. 



(ii) Condition 1 of the theorem ensures that the moment cone M,„ of (jA.ip is m-dimensional, i.e. not 
contained in some lower-dimensional hyperplane. 



Proof of Theorem \A.1[ Define the enlarged moment cone 

L ( /(x) 



M=ih = {h,...,bm+iV ■■ f f )Mdx), A*eM 
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and let M denote its closure. Then M is a closed convex cone and, moreover, the vector (&i, . . . , bm, V{D)) 
lies on its boundary. There exists a supporting hyperplane to M through this vector, specified by real 
constants {zi}'"'*'^ not all zero, such that 

rn 

(A.2) ^ + z„,+il/p) = 

i=l 

m+1 

(A.3) and ^ zA > 0, Vb e M. 

i=l 

In particular, on considering Dirac measures, one has 

m 

(A.4) ^'«*(^) + > 0, Vx e K. 

i=l 

We now show that z,„+i > 0. Indeed, for any 6 > 0, the vector {bi, . . . ,bm,V{D) — S) lies in the 
half-space complimentary to (jA.3|) . Therefore 

rn 

J2z^b^ + Z,n+l{V{D)-S)<0, V(5 > 0. 
1=1 

This clearly implies Zm+i > 0. However, z,„+i = is not possible, since this would contradict the 
assumption that b lies in the interior Mm, which is m dimensional. Thus, it must be that z,n+i > 0. 
Then from (IA.4I). it follows that 



/(x) > ^y,aj(x), 

1=1 

where yi := —Zi/zm+i for i = 1, . . . , m. Then (jA.2p becomes 

m 
i=l 

and this completes the proof. □ 

B Pathwise stochastic calculus 



This section describes a non-probabilistic approach to stochastic calculus, due to iFoUmen (Il98l[ ). that 
will enable us to define the continuous-time limit of the finite sums defining the realized variance, 
without having to assume that the realized price function < i— > 5*4 is a sample function of a semimartingale 
defined on some probability space. 

For T > let T = [0, T] with Borel sets B-y- A partition is a finite, ordered sequence of times 
■n = {Q = to < ti < ■ ■ ■ < tk = T} with mesh size m{'K) = maxi<j<fc(<j — ij-i). We fix a nested sequence 
of partitions {7r„, n = 1, 2, . . .} such that lim„_>.oo m(7r„) = 0. All statements below relate to this specific 
sequence. An obvious choice would be the set of dyadic partitions t" = jT/2" , j = 0, . . . , 2", n 1, 2, 

Definition B.l A continuous function X : T — !• K has the quadratic variation property if the sequence 
of measures on (T, 6t) 

tin^Y. {X{t, + l)-X{t,)f5t, 
tj £Tr„ 

(where St denotes the Dirac measure at t) converges weakly to a measure fi, possibly along some sub- 
sequence. The distribution function of fi is denoted {X)t, and we denote by Q the set of continuous 
functions having the quadratic variation property. 



Theorem B.2 ( (|Follmerl Il98lh ) For X e Q and f e C^, 



(B.l) f{Xt)-f{X,)= f'{Xs)dX, + - f"{X,)d{X) 



t 1 



2 Jo 
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Proof. (Sketch). The proof proceeds by writing the expansion 

(B.2) f{Xt) - fiXo) = Y.f\Xt^){Xt^^, - Xt^) + ^Y.f''(^t.)(^t.^^ - ^h)' + T.R^^^^ 

J j j 

where R{a,b) < (p{\b — a\){b — a)^ with (j) an increasing function, 0(c) — > as c — ^ 0. As m{Tr) -> 0, 
the third term on the right of (|B.2p converges to and the second term converges to the second term in 
(|B.1[) . This shows that the first term converges, so that essentiaUy the 'pathwise stochastic integral' is 
defined by (|B.1[) for arbitrary functions /'. □ 

An important remark is that Xt, being continuous, achieves its minimum and maximum X^ , X * in 
[0,T], so only the values of / in the compact interval (depending on the path X) are relevant 

in (IEtI) . 

For the applications in Section 0] we need an Tto formula' valid when /" is merely locally integrable, 
rather than continuous as it is in (jB.ip . In the theory of continuous semin iartingales, such an extensio n 
proceeds via local time and the Tanaka formula — see Theorem V L1.2 of Rogers an d Will iamsl ( 2000l ). 
Here we present a pathwise version, following the diploma thesis lUniversitat Bonn (,1980i ). For a 
function / we have f{b) — f{a) + J^^ f'(i/)dy, and applying the same formula to /' we obtain 



f{b)^f{a) = J (^/'(«)+y f"{u)dujdy = f'{a){b^a)+J{b-u)f"{u)du 

/OO 
i[aAb.aVb]iu)\b - u\f"{u)du. 
-OO 

Hence for any partition tt — {tj} of [0,t] we have the identity 

(B.3) f{Xt) - f{Xo) = nXt^)iXt,+. -Xt,)+ (l[x™,A-^=<](«)l^t.+. - «l) f'idu), 

where Xf'" = miii{Xt^, Xt^^,} and Xf'^'' = max{Xt^ , X*^^ J. We define 
(B.4) = 2^ l[;,™^;,_](u)|Xt,+, - u\, 

3 

and note that (u) = for u <^ {X^,X*). 

Definition B.3 Let C be the set of continuous paths on [0,T] such that the discrete pathwise local time 
Li" (u) converges weakly in L^(du) to a limit Lt{-), for each t € [0,T]. 

Proposition B.4 £ C Q. For X ^ C and t G [0, T] we have the occupation density formula 
(B.5) / Lt{u)du = [ lAiX,)d{X)s, A e 



Proof. Suppose X £ C, and let 7r„ be a sequence of partitions of [0,T] with m(7r„) — > 0. From (jB.3p we 
have, for f £ and t e [0,r], 

(B.6) fiXiJ ~ fiXo) - ^/'(XtJ(Xi^^, - Xt^) = - / Ll-iu)nu)du, 

where tn is the nearest partition point in 7r„ to t. As n — > cxd, the first term on the left of ()B.6p 
converges to f{Xt) and, since /" £ L^([ArH., X*], du), X G £ implies that the right-hand side converges 
to (1/2) / Lt{u)f"{u)du. Hence the 'integral' term on the left of (jB.6p also converges to, say, Lt. If we 
take f{x) ~ then the left-hand side of (|B.6p is equal to 

j:tj + l<t„ j:*3 + l<*n 

showing that (a) AT G Q and (b) Lt = /q f'{Xs)dXs, the Fdllmer integral of Theorem lB.2l It now follows 
from (lR6l) that 



(B.7) / Lt{u)f'{u)du^ /*/"(A,)d(A}„ 
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i.e., for any continuous function g we have 

Lt{u)giu)du= I g{X,)d{X), 



Approximating the indicator function 1^ by continuous functions and using the monotone convergence 
theorem, we obtain (IB.5I). □ 



For the next result, let W2 be the set of functions / in C^(]R) such that /' is weakly differentiable 
with derivative /" in Lf^^(R). 



Theorem B.5 If X £ C the Fdllmer integral extends in such a way that the pathwise Ito formula ()B.ip 
is valid for f G W2 . 

Proof. Let (phe a mollifier function, a non- negative C°° function on M such that (p{x) = for |a;| > 1 and 
/ 4>{x)dx = 1, define 0,i(a;) = n(f>(nx) and let fn{x) = / f{x — y)(j)niy)dy. Then /„, converge pointwise 
to /, /' and /" f" weakly in L^. From Proposition IB. 41 we know that X G Q and as /„ G we have, 
using (jB.ip and (|B.7() . 

fn{Xt) - fniXo) = / ./;(X,)dX, + - / Lt{u)f::{u)du. 
Jo ^ J-00 

As n 00, the left-hand side converges to f{Xt) — f{Xo) and the second term on the right converges to 
3 I-oo Lt{u)f"{u)du, and we can define 

f:AXs)AX, = hm \fn{Xt) - /„(Xo) - \ [ Lt{ii)f^{u)du 

It now follows from Proposition IB .41 that the Ito formula (jB.ip holds for this extended integral. □ 
Wc need one further result. 

Proposition B.6 Let X £ C and let f : M. ^ M. be a monotone function. Then Y — f{X) G C and 
the pathwise local times are related by 

(B.8) L^{u)^\f\f-\u))\Lf{f-\u)). 

Proof. We assume / is increasing — the argument is the same if it is decreasing — and denote v = f^^{u), 
so u — f{v). For the partition 7r„ we have for fixed t G T, from (|B.4p . 

if'""!") = IIl{/(A-™)</(.)</(A— 

3 

= E itc)<.<x»»)}I/'(-)(^*?,, - V)) + InoiXtj,, ~ vr\ 

(B.9) = /'(«)Elw-)<-<^-)} 



for ^ between 1; and Xt"^_^. Noting that is bounded for £^,v G [X^,X*] and that 

lim„^oo maxj \X(tj-_^_^) — X{tj-)\ = 0, we easily conclude that ii X £ C then the expression at (|B.9p 

converges in L^(dij) to f'{v)L^{v), so that Y E C with local time given by (|B.8[) . □ 

We note that in the above result it suffices to assume that / is defined on [X*, X*]. This allows us to 
apply the Proposition for / = log and and stock price trajectories in Section |4l Indeed, from (jB.8|) we 
obtain the elegant formula 

Lt°«^(«) = e-"Lf(e«). 

Finally, we build the connection between the pathwise calculus and the classical Ito (stochastic) 
calculus. 

Theorem B.7 Let {Xt,t G [0,T]) be a continuous semimartingale on some complete probability space 
(O, J^, P). Then there is a set N E T such that ¥N = and for lu ^ N the path t t— > X{t,Lu) belongs to 
C ( and hence to Q ), and Lt (u) defined in Definition \B.3\ coincides with the semimartingale local time of 
Xt at u. 
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Proof. First, Xt ^ Q a.s. Indeed, Theorem IV. 1.3 of Revuz and Yoil ( 1994 ) asserts that discrete approx- 



imations to the quadratic variation always converge in probabihty, so a sub-se quence c o nverg es almost 
surely, showing that £ Q in accordance with Definition IB. II It is shown by iFollmeil (|l98l[ i that the 



Ito integral and the pathwise integral defined by (jB.l[) coincide almost surely. Every semimartingalc S, 



3* 



has an associated local time which satisfies the occupation density formula (jB.Sp . It remains to show that 
with probability 1 th e discrete approximation s defined by (jB.4|) converge in L^(du). This is proved 
in the diploma thesis Universitat BonnI ( 1980l ). by detailed estimates which we cannot include here. □ 
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